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Topological crystalline superconductivity in the locally non-centrosymmetric multilayer supercon¬ 
ductors (SCs) is proposed. We study the odd-parity pair-density wave (PDW) state induced by the 
spin-singlet pairing interaction through the spin-orbit coupling. It is shown that the PDW state 
is a topological crystalline SC protected by a mirror symmetry, although it is topologically trivial 
according to the classification based on the standard topological periodic table. The topological 
property of the mirror subsectors is intuitively explained by adiabatically changing the BdG Hamil¬ 
tonian. A subsector of the bilayer PDW state reduces to the two-dimensional non-centrosymmetric 
SC, while a subsector of trilayer PDW state is topologically equivalent to the spinless p -wave SC. 

Chiral Majorana edge modes in trilayers can be realized without Cooper pairs in the spin-triplet 
channel and chemical potential tuning. 

PACS numbers: 74.20.Rp, 74.45.+C, 74.78.Fk 


Topologically nontrivial phases of superconductors 
(SCs) have evolved into one of the major research top¬ 
ics of modern condensed matter physics recentlyi^— . 
A characteristic feature of topological SCs is the fully 
gapped bulk spectrum accompanied by topologically pro¬ 
tected gapless edge states. Many of the topological su¬ 
perconducting states are realized in odd-parity SCs, and 
one of the most extensively studied examples is the chi¬ 
ral p x ± ipy-wave SC—However, only few materials are 
considered as possible hosts of odd-parity superconduc¬ 
tivity, because the conditions for spin-triplet pairing are 
quite unfavorable in most cases. So far, S^RuO^ii and 
some uranium-based heavy fermion compound s 18 : 19 show 
strong evidence for the spin-triplet odd-parity supercon¬ 
ductivity, but unfortunately their superconducting gap 
might have nodes on the Fermi surface. Recently, odd- 
parity topological superconductivity in a doped topologi¬ 
cal insulator Cu^B^Sea has been propose d 13 : 16 , however, 
experimental results are under debat e 20 : 21 . 

In a recent study we showed that odd-parity super¬ 
conductivity occurs naturally in multilayer systems with 
layer-dependent spin-orbit coupling arising from the lo¬ 
cal lack of inversion symmetry— We will consider here 
such locally non-centrosymmetric systems composed of 
the blocks of superconducting layers, e.g. trilayer sys¬ 
tems as depicted in Fig. [T] Here the layer-dependent 
Rashba spin-orbit coupling is responsible for unusual 
electronic and superconducting properties— The cou¬ 
pling constant of the Rashba spin-orbit coupling shows 
the layer-dependence, (ay,CC 2 ,<^ 3 ) = (a,0, — a), ensured 
by the global inversion symmetry. We have shown that 
in such a system an odd-parity superconducting state 
can be stabilized by a magnetic field, even if the zero- 
field phase is the even-parity state (see Fig. [U -■ To 
be precise, the order parameter in the spin-singlet chan¬ 
nel changes sign between the outer-most layers in the 
field-induced superconducting state (see Fig. [[])• Con¬ 
sidering the spatially modulating order parameter in 
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FIG. 1. (Color online) Schematic figure of trilayer system. 
The filled (open) circles represent the 2D superconducting 
(normal spacer) layers. The dashed line denotes the mir¬ 
ror plane. Attached fists provide information on the layer- 
dependent Rashba spin-orbit coupling and the order parame¬ 
ters in the BCS and PDW states. 


the trilayer, we call it the “pair-density wave (PDW) 
state”—. Multilayer structures of this kind are not only 
theoretical constructs, but have indeed been produced 
recently, for example, in the artificially grown super¬ 
lattices CeCoIns/YbCoIn^F— and in transition-metal- 
oxide interfaces— The PDW state is stabilized when 
the three conditions, (a) Pauli-limited SC, (b) quasi-two- 
dimensional structure, and (c) large spin-orbit coupling, 
are satisfied. These conditions are naturally satisfied in 
the heavy fermion superlattice CeCoIns/YbCoIng^iT— . 
Furthermore, the recent technology enabled the artificial 
tuning of the superlattice structure— Thus, we may 
expect that the PDW state can be stabilized in a su¬ 
perlattice CeCoIns/YbCoIns, although no experimental 
evidence has been reported so far. In view of the experi¬ 
mental and theoretical status, the discussion of topolog¬ 
ical features of the PDW state is well motivated. 

Topological aspects of bilayer PDW state in the ab- 
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sence of a magnetic field have been investigated by Nako- 
sai et al£L. They showed that the bilayer PDW state is a 
topological state protected by a Z 2 invariant when (and 
only when) the Fermi level lies in the hybridization gap 
between the bonding and anti-bonding bands. The field- 
induced PDW phase in the multilayer system has not 
been investigated in this respect so far. 

First, we consider the topological properties of the 
PDW state on the basis of the so-called topological peri¬ 
odic tabled. When time-reversal symmetry is broken by 
a magnetic field, the symmetry class of the state is D. 
The two-dimensional (2D) system in the class D is char¬ 
acterized by an integer topological number, the Chern 
numbe r 30 ! 31 . However, the Chern number must be zero 
in the time-reversal invariant system, and the magnetic 
field does not change the Chern number without clos¬ 
ing the gap. According to the numerical analysis of 
the Bogoliubov-de-Gennes (BdG) equation, the magnetic 
field does not close the gap in the PDW stated. Thus, 
this shows that the field-induced PDW state is topolog¬ 
ically trivial in terms of the classification based on the 
topological periodic table. 

On the other hand, recent developments in the clas¬ 
sification scheme of topological phases shed new light 
on topological phases protected by the crystal sym- 
metrj*22r— . The ’’topological crystalline SCs” have 
been classified relying on the mirror, inversion, rotation 
and magnetic point group symmetr y 35 ! 37 ' 39 . The spin- 
triplet superconducting/superfluid states in S^RuOzt^li, 
UPtg—, and 3 He22 have been discussed from this point 
of view. In this letter, we will show that the spin-singlet 
PDW state in trilayers is generally a topological crys¬ 
talline SC protected by the mirror symmetry. This is, 
to our knowledge, the first proposal for the topological 
crystalline SC without requiring the pairing interaction 
in the spin-triplet channel. 

We consider the mean-field BdG Hamiltonian for the 
2D multilayer SC, 

fi= ^2 [£(fc)(7 0 + Ot m g{k) ■ (T - HBH(J z ] ss 'C' ksrn C ks 'm 

k,s,s' ,m 

E C ksm C ksm' 

k,s,(m,m') 

k,s,s' ,m 

where k, s, and to (= 1 are indices of mo¬ 

mentum, spin, and layer, respectively. We assume the 
simple dispersions £(fc) = ~2t(cosk x + cos k v ) — /i and 
g(k) = (— sin k v , sin k x , 0). The latter describes the 
Rashba spin-orbit coupling, whereby the coupling con¬ 
stant a m is layer-dependent. Nearest-neighbor layers 
are coupled by the hopping matrix element t±. We fo¬ 
cus on the intra-layer Cooper pairing which is relevant 
for 2D SCs, as realized in CeCoIn 5 /YbCoIn 5 superlat- 
tices^r— and <5-doped SrTiO^, although an inter-layer 
Cooper pairing has been considered for Cu^BiaSe a 13 ! 16 . 


The layer-dependent order parameter can then be pa¬ 
rameterized by A m (k) = [i/j m (k) + d m (k ) • a]ia y , where 
'tp m (k) and d m {k) represent the spin-singlet and spin- 
triplet components of order parameters on the layer to, 
respectively. For simplicity, we assume the S+p -wave 
pairing state, in which the dominant s-wave order pa¬ 
rameter ip m (k) = if m is mixed with the spin-triplet p- 
wave component through spin-orbit coupling and pair¬ 
ing interaction. The latter has the structure d m (k) = 
G m (— sin k y , sin k x , 0) + ib m ( sin k x , sin k y , 0), obtained by 
solving the BdG equation^. In the following we ana¬ 
lyze the two competing solutions of the BdG equation: 
(1) the ”BCS state” with f} m {k) = ipM+i-m{k) and 
d m (k) = —dM+i-m{k) and (2) the ”PDW state” where 
V’m(fe) = - ipM+i-m(k ) and d m (k ) = d M +i-m(k). We 
now assume a pairing mechanism favoring spin-singlet 
pairing, as often given by electron-phonon coupling or 
antiferromagnetic spin fluctuation. Thus, the BCS state 
is stabilized by the inter-layer Josephson coupling at zero 
magnetic field. However, the PDW state is stabilized by 
spin-orbit coupling in the high magnetic field region at 
sufficiently low temperatures^. 

Now we define the topological invariant of multilayer 
SCs protected by the mirror symmetry, by means of the 
mirror Chern number. The BdG Hamiltonian is repre¬ 
sented as ,'H = \ ]T) fc 'Fj c H(fe)'I , fc with use of Nambu op¬ 
erators 'kj, = ( c i S m’C-ksm ) in 4 x M dimension. The 
mirror symmetry with respect to the central a;?/-plane is 
obeyed, 

M%U{k)Mtt =U{k). (2) 

A4jjr y is the mirror reflection operator in the particle-hole 
space (see Appendix A). We introduce M xy for the BCS 
state and Mf y for the PDW state, respectively. Equa¬ 
tion m guarantees that the BdG Hamiltonian can be 
block-diagonalized in the eigenbasis of M xy . Thus, the 
system is divided into the two subsectors corresponding 
to the block Hamiltonian "H^fc) with A = ±i as eigen¬ 
values of We now define the mirror Chern number 

z/(A), as the Chern number of the subsector Hamilto- 
nia n 33 ! 36 (see Appendix B). The topological protection 
of the mirror Chern number is guaranteed in some topo¬ 
logical classes characterized by the symmetries of sub¬ 
sector Hamiltonian 'H±(k)r. Important here are the 
time-reversal, particle-hole and chiral symmetry (see Ap¬ 
pendix C). 

For illustration we first discuss the bilayer system. We 
obtain the subsector Hamiltonian for the X = i sector as 

± / W(k) +t ± a z ±i[ip - d(k) ■ cr]a v A 

A= * \ ^ia y [il)* — d*{k) ■ <x] — Ti T (—k) ± t±a z J ’ 

( 3 ) 

where "H'(fc) = ^(fc)cro — pbHct z — ag(k) ■ <j. The sub¬ 
sector Hamiltonian for A = — i is obtained by changing 
the sign of t±, as t± —t±. For the BCS state, al¬ 
though the particle-hole symmetry in the original BdG 
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Hamiltonian is conserved, we cannot rely on this sym¬ 
metry in the subsector Hamiltonian unless the special 
condition namely t± = 0, is satis¬ 

fied (demonstrated in the supplementary material). On 
the other hand, the chiral symmetry is conserved in this 
subsector at H = 0. Therefore, in the absence of a mag¬ 
netic field the symmetry class is AIII which is topologi¬ 
cally trivial in 2DA If the chiral symmetry is broken by 
a magnetic field, both subsectors belong to the class A, 
which is characterized by an integer topological invari¬ 
ant^. However, both subsectors are topologically trivial, 
i/(A) = 0, or the gap is closed under the realistic condi¬ 
tion, \ip\ < t±. 

For the odd-parity PDW state, time-reversal sym¬ 
metry in the subsector Hamiltonian is ill-defined for 
t± ^ 0, while the particle-hole symmetry is con¬ 
served. Thus, the subsector belongs to the symme¬ 
try class D unless ( t±,H ) = (0,0). Interestingly, 
each subsector is equivalent to the BdG Hamiltonian 
of a 2D non-centrosymmetric superconductor (NCSC)^ 
with the fictitious magnetic field H ± t±, whose 
topological property has already been clarified^ -10 ’ 12 . 
The dominantly spin-singlet pairing state |d(fe)| < 

\ip\ can be topologically nontrivial, when the effec¬ 
tive magnetic field /ibH ± tj_ satisfies the condition 
v/ (4t + /x) 2 + \ip\ 2 < \fi B H±t±\ < ^(4t - n) 2 + \i/j\ 2 , 

IvV + IV'I 2 < < yj (4t - /i) 2 + |V>| 2 ] for 

[i < — 2t [— 2t < /i < 0J— Although great effort has 
been devoted to the realization of this condition in semi¬ 
conductor devices^, this condition needs fine tuning of 
the chemical potential and is rather unrealistic in metals. 

For H = 0, this condition is indeed equivalent to 
the criterion for a Z 2 topological SC without relying on 
the mirror symmetry— This means that the nontrivial 


Z 2 topological number in the original BdG Hamiltonian 
(class Dill) is obtained by the mirror Chern number of 
the subsectors (class D). This is analogous to the fact 
that some Z 2 topological insulators are characterized by 
the spin Chern number—. Our analysis sheds light on 
the analogy between the 2D NCSC and the Z 2 nontrivial 
bilayer SC, the former being equivalent to a mirror sub¬ 
sector of the latter. The interlayer coupling t± plays the 
same role as the magnetic field in the former. Although 
the Z 2 number of the original BdG Hamiltonian is not 
a topological invariant in the presence of the magnetic 
field, the mirror Chern number is topologically protected. 
Therefore, the mirror Chern number is useful to indicate 
the topological property of field-induced superconducting 
states. 

We now turn to the trilayer system to show the most 
important results of this paper. We consider the trilayer 
structure conserving the mirror symmetry (see Fig. [T|), 
and adopt the layer-dependent Rashba spin-orbit cou¬ 
pling ( 0 : 1 , 012 , 0 : 3 ) = (o, 0, —o). The layer-dependent or¬ 
der parameters are shown in Fig. [0 Using the mirror 
operator with respect to the central xy-plane, the BdG 
Hamiltonian is again block-diagonalized into the mirror 
subsectors. We show the subsector Hamiltonian for the 
BCS state in Appendix B. The subsector belongs to the 
class A for H ^ 0 and to the class AM for AT = 0, if 
t± ^ 0. We confirmed that the mirror Chern number is 
zero or the gap is closed as in bilayers. Thus, topolog¬ 
ical superconductivity is not realized in the BCS state. 
Indeed, Fig. EJa) shows no zero energy Majorana mode, 
indicating the topologically trivial property. 

In contrast, the PDW phase represents a topological 
crystalline superconducting state. We obtain the subsec¬ 
tor Hamiltonian 


/ 

£t(*0 

ak+ 

V^L 

0 

^out— (^) 

b'oiit 

\ 


ak- 

a(*o 

0 

0 

V^out 

dout+ (k) 



t± 

0 

£t(*0 

din— (&) 

0 

0 



0 

0 

-d* n .(k) 

-ftW 

H 
. ■ + ° 

*> 

1 

0 



~d* 0 ut-(*0 

4>out 

0 

—\ 


ak- 


V 

-V’out 

^out+(^) 

0 

0 

ak+ 

-a(fc) 

/ 


for A = i. We denote £ s (fc) = £(fc) - (<j z ) ss hbH, 
k ± = sin k y ± i sin k x , and d 0 ut(in)±( fc ) = 4 ut(in)( fc ) ± 
*^out(ir ^ (k). The subsector Hamiltonian for A = —i is 
shown in the supplemental material. Both subsectors 
belong to the symmetry class D independent of the mag¬ 
netic field, if t _l ^ 0. Therefore, the mirror Chern num¬ 
ber is a topological invariant. We obtain a nontrivial 
mirror Chern number zz(A = ±z) = =pl, almost indepen¬ 
dent of the parameters. In contrast to the bilayer PDW 
state, this topologically nontrivial superconducting state 
is realized without having to rely on a special choice of 


parameters. Because the mirror Chern number is odd, 
the trilayer PDW state is also a Z 2 topological super¬ 
conducting state at H = 0, although the magnetic field 
is required for the thermodynamic stability of the PDW 
state— 

An intuitive understanding of our result can be ob¬ 
tained by adiabatically deforming the subsector Hamil¬ 
tonian Hyik). The interlayer coupling t± is decreased 
to zero without closing the gap as long as the spin-triplet 
component di n (k) is finite. The topology does not change 
through this adiabatic deforming. Then, the finite mir¬ 
ror Chern number originates from the decoupled 2x2 
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matrix in the center of 6 x 6 matrix of Eq. which 
denotes a spinless chiral p- wave SC. It has been shown 
that the spinless chiral p-wave SC is topologically non- 
triviali and the Chern number is ±1 (see Appendix D). 
Indeed, we obtained the nontrivial mirror Chern number 
zz(±*) = =Fl, which is identified as the Chern number 
originates from the decoupled 2 x 2 matrix in the limit 
t± —> 0. Now it became apparent that no fine tuning 
of the chemical potential is needed. The other 4x4 ma¬ 
trix decoupled in the subsector Hamiltonian describes the 
2D Rashba-type NCSC which has been proposed to be a 
topological s-wave SC^ d°i 12 . However, we do not assume 
a fine tuning of the chemical potential which is required 
in their proposals. 

We emphasize that the Cooper pairing in the p- wave 
channel d ln / out (k) is not needed for the topological crys¬ 
talline superconductivity, although it played an impor¬ 
tant role in the above intuitive explanation. This is un¬ 
derstood from the fact that d in i ont (k) is decreased to 
zero without closing the gap when the interlayer hopping 
f_L is finite— Thus, the topology is equivalent between 
the Hamiltonian for t± = 0 and d ln / ont (k) ^ 0 (as in 
the above intuitive explanation) and that for t± ^ 0 and 
din/out(fc) = 0 ( as we consider here). This means that the 
topological crystalline superconductivity is realized with¬ 
out any attractive interaction in the spin-triplet channel. 
Once the PDW state is stabilized in the trilayer system, 
it is a topological crystalline SC. 

In order to verify the bulk-edge correspondence, we 
show the presence of edge states in the trilayer SCs. Fig¬ 
ures m a) and (b) show the energy spectra of BCS state 
and PDW state, respectively, for a ribbon-shaped system 
with open boundaries along x-axis and translational in¬ 
variance along y-direction. Consistent with the vanishing 
mirror Chern number, no subgap edge state appear in the 
BCS state. In contrast, we find two chiral Majorana edge 
modes in the PDW state. One comes from the A = i sub¬ 
sector (solid lines) and the other comes from the A = — i 
subsector (dashed lines). These modes are not Kramers 
pairs, because the time-reversal symmetry is broken by 
the magnetic field. We confirmed that the presence of 
these Majorana modes is robust against the change of 
parameters, such as variations of V’in/out; ain/out, &m/out> 
t_ l, a, and p. 

In Figs. [21(c) and (d), we show the spatial profiles of the 
zero-energy Majorana modes localized around the edge. 
Large probability density on the inner layer \<fi S 2(x)\ = 
\(x,s2\E = 0)| is also shown. This means that the Ma¬ 
jorana state mainly originates from the inner layer, as 
expected from the intuitive explanation discussed above. 

In this letter we have focused on the 2D multilayer 
SCs, but the topologically nontrivial properties also ap¬ 
pear in the three-dimensional (3D) system. When we 
take into account an inter-multilayer coupling through 
normal spacer layers (see Fig. [T]) and consider the 3D 
Brillouin zone, the BdG Hamiltonian conserves the mir¬ 
ror reflection symmetry as A4± y 'H(k x ,k y ,k z )Ai xy = 
7-L(k x ,k y ,—k z ). Thus, the mirror symmetry defined by 



FIG. 2. (Color online) Energy spectra of (a) the BCS state 
and (b) the PDW state with open boundaries at x = 1 and 
x = 200. The solid and dashed lines in (b) show the Majorana 
edge modes in A = i and A = — i subsectors, respectively. 
Thick (green) lines show the edge states near the boundary 
x = 1, while thin (red) lines show the edge states near x = 
200. We take t = 1, p = —2, pbH = 0.3, a = 0.3, t± = 0.1, 
Ipout — 'dirI — 0.5, (Tout — din = 0.05, and &out — 5in — 
0.1. (c) and (d) illustrate the wave function of Majorana 

modes localized around x = 1. Amplitude of spin- and layer- 
resolved wave function, (p S m(x ) = (x,sm\E = 0), is shown. 
The Majorana state resides dominantly on the center layer 
(m = 2 ) with up spin for the subsector A = i (c) and with 
down spin for A = — i (d). 


Eq. © is satisfied in the mirror invariant planes, k z = 0 
and 7 r. We can define the mirror Chern number in these 
2D mirror invariant planes, and we indeed obtain a non¬ 
trivial mirror Chern number at both k z = 0 and 7 r for a 
small inter-multilayer coupling. We confirmed that Ma¬ 
jorana cones appear on [ 100 ] and [ 010 ] surfaces where the 
mirror symmetry is conserved. 

Analyzing topological properties of multilayer SCs we 
found that the PDW state is a topological crystalline su¬ 
perconducting phase protected by the mirror symmetry. 
We stress that a purely s-wave PDW state in trilayers 
can be a topological SC accompanied by the Majorana 
fermion on its edge without the tuning of chemical po¬ 
tential, which is necessary in the bilayer PDW stated 
and the ID and 2D NCSCSi— This finding significantly 
expands the possibility of realizing the topological SC be¬ 
cause most SCs have a s-wave symmetry. It is straight¬ 
forward to extend our analysis to more than three layers 
and we find that the PDW state is a topological SC in¬ 
dependent of parameters, if the number of layers is odd. 
Thus, the design of the topological crystalline SC is fea¬ 
sible for artificially grown multilayers using the available 
technology^—. The superlattice CeCoIns/YbCoIn^ 
is considered to be a D+p -wave SC, and will be sim¬ 
ilarly a topological crystalline SC as will be discussed 
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Appendix A: MIRROR OPERATOR FOR 
MULTILAYER SUPERCONDUCTOR 


In this section, we introduce the mirror operator for 
multilayer SCs. For the clarity of discussions, we de¬ 
scribe the BdG Hamiltonian with use of the normal state 
Hamiltonian Ho{k) and the gap function A(fc) as, 


H{k) = 


'Hoik) A (fc) A 

At(fc) -nli-k)J ■ 


(Al) 


Because the mirror symmetry with respect to the xy- 
plane is evidently conserved in the multilayer systems, 
the normal state Hamiltonian 'Hoik) is invariant for the 
mirror reflection operator M xy as, 


M xy Uo{k)Ml y =Uo{k). (A2) 

The mirror reflection operator transforms the momentum 
k as, 


k — ikxi k y , k z ) ^ ikxi k y , k z ), (A3) 

and the spin as, 

S — ($xi Sy-i %z) t ( S x , Sy, S^,), (A4) 

respectively. As we focus on the 2D system, the momen¬ 
tum is invariant under the mirror reflection, while the 
spin is subject to the 7r-rotation around the z axis. In 
addition to these fundamental degrees of freedom, the 
multilayer systems also have the index for layer to. By 
taking into account the reversal of layers, the mirror re¬ 
flection operator for the normal part is given by 


A d xy — 


ia z 0 
0 ia z 



0 ia z 
ia z 0 


, (A5) 

I 


for bilayers, while it is given by 

/ ic z 0 0 \ / 0 0 a 0 \ 

M xy = 0 ia z 0 x 0 (T 0 0 

\ 0 0 ia z J \ (T 0 0 0 ) 

( o 0 ia z \ 

= 0 ia z 0 , (A6) 

\ia z 0 0 J 


for trilayers. It is straightforward to obtain the mirror 
reflection operator for more than three layers. We con¬ 
firmed that the normal state Hamiltonian Ho(k) is in¬ 
variant for these operators. 

When the gap function has a well-defined mirror-parity 
as M xy Aik)M xy = ±A(fe), the BdG Hamiltonian H{k) 
is invariant under the mirror reflection as, 

M%H{k)M% = H{k), (A7) 

where the mirror operator in the particle-hole space M xy 
is introduced as, 


M 


± 

xy 


M X y 0 \ 

0 ±M* xy J ■ 


(A8) 


The sign + (—) is adopted in the even-parity (odd-parity) 
superconducting state. Thus, we adopt M xy in the BCS 
state while M~ in the PDW state. 


Appendix B: MIRROR CHERN NUMBER 

Because the BdG Hamiltonian H(k) and the mirror op¬ 
erator M xy are commutative, we can block-diagonalize 
the BdG Hamiltonian using the eigenbasis of mirror op¬ 
erator, 

^ w± , = (wiA> w£ ° ((fc) ). (Bl) 

The unitary matrix U ± is obtained by the eigenbasis of 
M xyl and the subsector Hamiltonian is characterized by 
the eigenvalues, A = Hi. Examples of subsector Hamil¬ 
tonian for the bilayer system and the trilayer PDW state 
for A = i are given in Eqs. (3) and (4). The subsector 
Hamiltonian for the BCS state in trilayers is given by 


/ £t(*0 

ak+ 

V2tj_ 

0 

d 0 ut-ik) 

V^out 

\ 


a(*o 

0 

0 

d’our 

^out+ (^) 



0 

et(fc) 

Ipin 

0 
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0 

0 

Cn 

-a(fc) 

0 

H 

, -to 

i 



-Cut 
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-£t(*0 



V V'out 

— ^out +(k) 
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— \/2t± 

ak+ 

-a(fc) 
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(B2) 


and 
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o = 


( £t(*0 

ak + 

0 

0 

dout- (^) 

Pont 

\ 

ak- 

Ct(fc) 

V2t ± 

0 

Pout 

dout-f- (k) 


0 

a/2 t± 

a(*o 

-V’in 

0 

0 


0 

0 

-c* 

Ct(*0 

—a/2 t± 

0 


d* ou t_(fc) 

“Cut 

0 

— a/2 t± 

-w 

ak- 


V Cut 

— d* ut +(^) 

0 

0 

ak+ 

-&(*) 

/ 


(B3) 


r 


while we obtain the A = 
the PDW state as 


—i subsector Hamiltonian for 


Kx = _ t (k) = 


/ 

£t(fc) 

ak+ 

0 

0 

dout— (^) 

Pou t ^ 


ak- 

a(fc) 

V2t ± 
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Pout 

dout+ (^) 
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A/2t_L 

£l(*0 

din-f- (^) 

0 
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dfn+W 
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— a/2*_l 


-d* ut _(fc) 

Cut 
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Ct(fc) 

ak _ 
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-Cut 

dout+( fc ) 
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aP+ 

-&(*0 / 


(B4) 


The mirror Chern number za(A) is defined by 

KA) = ±ff dk x dk y e^d ki Af x (k), (B5) 

A t\( k ) =i ( U n,A( fc )l 5 fc« U n,A( fc ))’ ( B6 ) 

CC fc ) <0 

where E± x (k) and \u± x (k)) are the eigenenergy and 
eigenstate of the subsector Hamiltonian 7i x (k), namely, 

^(k)\u± x (k)) = £± A (fc)|< A (fc)>. (B7) 

For the numerical calculation of the mirror Chern num¬ 
ber, we adopt an efficient method developed in Ref. l46l . 


Appendix C: SYMMETRY CLASS OF THE 
SUBSECTOR HAMILTONIAN 


We examine the symmetry class of the subsector 
Hamiltonian T-L x (k). For this purpose, we first consider 
the time-reversal symmetry, particle-hole symmetry, and 
chiral symmetry in the original BdG Hamiltonian, which 
are defined as, 


TU{k)T f = H*{-k ), 

PH(fc)P t = -n T {-k), 
CH(k)C f = -P(fc), 


(Cl) 

(C2) 

(C3) 


respectively. For bilayers, the operators T, P, and C are 
given by 


T = diag(itJ y , ia y , ia y , ia y ), 
( 0 0 ff 0 0 

0 0 0 Co 

do 0 0 0 

V 0 (To 0 0 


P = 


(C4) 


(C5) 


I- 

and C = PT\ respectively. As a result of the uni¬ 
tary transformation with use of V , Eqs. (EH)-(PHI) are 
rewritten in terms of the subsector Hamiltonian. For the 
BCS state, we obtain 


i(j y 0 


o' JJ^)( ~o v 


0 


= n+l_ i (-k), 


(C6) 


U o)* 5 -. ( ‘>(° (CD 

) pm ( 7 ) = -pm- < C8 > 


0 — ia y 

0 


—la 


Equation (IC6I) indicates that the time-reversal symme¬ 
try in the subsector is ill-defined unless r H X:=i (k) = 
'H x= _ i {k) : even when the time-reversal symmetry is con¬ 
served in the original BdG Hamiltonian. Only when the 
special condition P A=i (fc) = P A= _ i (fc) is satisfied and 
Eq. El holds, we can rely on the time-reversal symme¬ 
try in the subsector Hamiltonian. Similarly, the condi¬ 
tion P A=i (fc) = % A= _j(fc) as well as Eq. (IC7I) have to be 
satisfied for the particle-hole symmetry in the subsector 
Hamiltonian. On the other hand, the chiral symmetry is 
well-defined in the subsector Hamiltonian as Eq. (IC8I) . 

In the same way, the time-reversal symmetry and the 
particle-hole symmetry in the PDW state are described 
with use of the subsector Hamiltonian as, 


0 


0 ia, 


'v 


v 





n- x (k) 


0 —ia, 




= -H7 1 (-k), (CIO) 
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t± = 0, H = 0 

tj ^ 0, H = 0 

t ± =0, H ^0 

otherwise 

BCS state 

Dill 

AIII 

D 

A 

PDW state 

Dill 

D 

D 

D 


TABLE I. Symmetry class of the subsector Hamiltonian T-L^{k) for the bilayer SC. 


respectively. According to Eq. m , the particle-hole 
symmetry is well-defined in the subsector Hamiltonian, 
and it is always conserved in the PDW state. Thus, we do 
not have to consider the chiral symmetry because it coin¬ 
cides with the time-reversal symmetry. The time-reversal 
symmetry is conserved in the subsector Hamiltonian only 
when the special condition T~L^ =i {k) = 'H A= _ i (fc) is sat¬ 


isfied and Eq. (IC10I) holds. In Table Q] we summarize 
the symmetry class of the subsector Hamiltonian for the 
bilayer BCS and PDW states. 

The symmetry class of the subsector Hamiltonian in 
the trilayer SC can be analyzed in the same way. For 
the BCS state, the time-reversal symmetry, particle-hole 
symmetry, and chiral symmetry are given by 


( -lOy 

0 

0 \ 

( 

lOy 

0 

° 


0 


O 

3i 

v + 

Jj. 


0 

cr z 


(Cll) 

V 0 

0 

iljy J 

\ 

0 

0 

ia y J 



/ 0 

0 


\ 

( o 

0 



0 

&X 

0 


0 

® X 

o \ =-ntl_ i {-k) 

(C12) 

\ cr 0 

0 

0 y 

/ 

\ cr 0 

0 

0 J 



/ 0 

0 

~ia y \ 

/ 

0 

0 

ia y \ 


0 

^i(Jy 

o n+ =±i {k) 


0 

±i<Ty 

o \=-n+ =±i (k), 

(CIS) 

\ ~ l(J v 

0 

0 / 

V 

ia y 

0 

0 J 




tx = 0, H arO 

t± 0, H = 0 

otherwise 

t± = 0, H = 0, di n = 0 

otherwise 

Symmetry class 

Dili 

AIII 

A Symmetry class 

Dill 

D 


TABLE II. Symmetry class of the subsector Hamiltonian TABLE III. Symmetry class of the subsector Hamiltonian 
'H'l(k) in the trilayer BCS state. in the trilayer PDW state. 


respectively. On the other hand, the time-reversal sym¬ 
metry and particle-hole symmetry in the PDW state are 
given by 


/ 

ICTy 

0 

° ^ 


ia y 

0 

° ^ 


0 

CTO 

0 H A=i 


0 


° 

V 

0 

0 

-lOy ) 

\ 

0 

0 

l(Jy ) 


(C14) 


l 0 

0 

era \ 

( 0 

0 

CTO \ 

0 

&x 

o \n^(k) 

0 

&x 

° 

V o-o 

0 

0 J 

\ ^0 

0 

0 ) 


(CIS) 


The particle-hole symmetry is always conserved in the 
subsector Hamiltonian of the PDW state. We summarize 
the symmetry class of the subsector Hamiltonian for the 
BCS and PDW states in Tables HD and Hill respectively. 
It is straightforward to elucidate the symmetry class 


of more than three layers. Independent of the number 
of layers, the subsector Hamiltonian belongs to the class 
D (class A) in the PDW state (BCS state), under the 
realistic conditions t± ^ 0 and ff / 0. 


Appendix D: TOPOLOGICAL NUMBER OF 
SPINLESS CHIRAL P-WAVE SC 

The topological number of spinless chiral p- wave SC 
part in Eqs. GD and (iBdl) is given by & 

= sgn^A(fc 0 )]sg n [det5 fei ??Aj(fc 0 )],(Dl) 

T 7 A (fc °)=0 

where £ A (fc) = £j(fe) and r]\(k) = —di n _(fc) for A = i, 
and £ A (fc) = £j,(fe) and ?y A (fc) = d in+ (k) for A = —i. 
We denoted r/\(k) = r]x,i{k) +i^\^{k). We find i/ p (A) 
to be non-zero if the odd number of zero-nodes of rj\(k) 











































[r]\(k°) = 0] are enclosed by the Fermi surface. Indeed, 
one zero node is enclosed by the Fermi surface for the 
simple dispersion adopted in this paper. Thus, we ob¬ 


tained the nontrivial Chern number v p (±i) = =p l, which 
is identified as the mirror Chern number of the subsector 
Hamiltonian Eqs. (j4j and (IB4I) in the limit t± —> 0. 
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